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ABSTRACT 

Spatial inhomogeneities in the spectral shape of the ultra-violet background 
(UVB) at the tail-end of He n reionisation are thought to be the primary cause of the 
large fluctuations observed in the Hen to Hi Lya forest optical depth ratio, THeii/THi, 
at z ~ 2 — 3. These spectral hardness fluctuations will also influence the ionisation 
balance of intergalactic metals; we extract realistic quasar absorption spectra from a 
large hydrodynamical simulation to examine their impact on intergalactic Siiv and 
C iv absorbers. Using a variety of toy UVB models, we find that while the predicted 
spatial inhomogeneities in spectral hardness have a significant impact on THeii/Tm, the 
longer mean free path for photons with frequencies above and below the He n ionisation 
edge means these fluctuations have less effect on the Siiv and C iv ionisation balance. 
Furthermore, UVB models which produce the largest fluctuations in specific inten- 
sity at the He n ionisation edge also have the softest ionising spectra, and thus result 
in photo-ionisation rates which are too low to produce significant fluctuations in the 
observed tsuv/^CIV- Instead, we find spatial variations in the IGM metallicity will 
dominate any scatter in tsuv/tcw- Our results suggest that observational evidence 
for homogeneity in the observed tsiiv/tciv distribution does not rule out the possi- 
bility of significant fluctuations in the UVB spectral shape at z ~ 2 — 3. On the other 
hand, the scatter in metallicity inferred from observations of intergalactic C iv and 
Si iv absorption at z ~ 2 — 3 using spatially uniform ionisation corrections is likely in- 
trinsic, and therefore provides a valuable constraint on intergalactic metal enrichment 
scenarios at these redshifts. 

Key words: methods: numerical - intergalactic medium - quasars: absorption lines. 



S-i 1 1 INTRODUCTION 

a ■ 

High resolution quasar absorption line spectroscopy has en- 
abled the statistical detection of intergalactic metals at cos- 
mic overdensities as low as A = p/{p) = 1 — 10 at z ~ 3. 
The inferred metallicities are around 10~ 2 — W~ 3 Zq, with 
C iv lines associated with a substantial fraction of Lya for- 
est lines with N m > 10 14 ' 5 cm -2 (Cowie et al. 1995; Ellison 
et al. 2000; Simcoe et al. 2004; D'Odorico et al. 2010). Ad- 
ditional metal lines, such as Siiv (Songaila & Cowie 1996; 
Aguirre et al. 2004) and Ovi (Schaye et al. 2000; Pieri & 
Haehnelt 2004; Aguirre et al. 2008) are also detected in the 
low density intergalactic medium (IGM). Observations indi- 
cate metallicity increases with density (Schaye et al. 2003; 
Aracil et al. 2004) and the distribution of metals in the IGM 
is inhomogeneous (Scannapieco et al. 2006; Pieri et al. 2006; 
Schaye et al. 2007; Fechner & Richter 2009; Martin et al. 
2010). 



A crucial ingredient in many of these studies are the ion- 
isation corrections which convert the observed ionic abun- 
dances to metallicities. These corrections depend on the den- 
sity and temperature of the IGM, as well as the intensity and 
spectral shape of the metagalactic ultra-violet background 
(UVB). The former may be modelled using cosmological hy- 
drodynamical simulations which follow the enrichment his- 
tory of the IGM (e.g. Theuns et al. 2002a; Cen et al. 2005; 
Oppenheimer & Dave 2006; Tescari et al. 2010). The lat- 
ter is usually calculated using detailed models for the UVB 
intensity and spectral shape (Haardt & Madau 1996, 2001; 
Fardal et al. 1998; Faucher-Giguere et al. 2009). Simulations 
of intergalactic metal enrichment typically assume the meta- 
galactic ionising radiation field is spatially uniform. This is 
expected to be a reasonable approximation for hydrogen ion- 
ising photons at z ~ 3, when the mean free path of ionising 
photons is much larger than the average separation between 
ionising sources (Meiksin & White 2004; Croft 2004). 
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However, the epoch of Hen reionisation is thought to 
complete around z ~ 3, when quasars are numerous enough 
to provide the hard photons (E > 4Ryd) needed to dou- 
bly ionise helium (Furlanetto & Oh 2008; McQuinn et al. 
2009). Indeed, observations of the Hi and Hen Lya forest 
indicate the UVB spectral shape fluctuates significantly on 
scales of 4 — 10 Mpc at these redshifts, producing a wide 
range of values for the Hen to Hi column density ratio, 

V = Nhcu/Nhi ~ 1 - 10 3 (Zheng et al. 2004; Fechner et al. 
2006; Shull et al. 2004, 2010). These fluctuations are thought 
to arise from the small number of quasars lying within the 
short (Hen ionising photon) mean free path expected at 
the tail-end of Hen reionisation (Fardal et al. 1998; Bolton 
et al. 2006; Furlanetto 2009a), although small scale radiative 
transfer effects (Maselli & Ferrara 2005; Tittley & Meiksin 
2007), collisionally ionised gas (Muzahid et al. 2010) or a 
significant number of thermally broadened lines (Fechner & 
Reimers 2007) may also produce column density ratios with 

V < io. 

These spatial inhomgeneities in the UVB spectral shape 
should also have an impact on the ionisation balance of the 
metals in the IGM. 1 Intriguingly, independent of the spec- 
tral shape of the (spatially uniform) UVB used to make ion- 
isation corrections, several studies find evidence for scatter 
in the IGM metallicity at fixed density (Rauch et al. 1997; 
Schaye et al. 2003; Simcoe et al. 2004). This scatter may 
be intrinsic, or it may instead be partially attributable to 
spatial fluctuations in the spectral shape of the UVB. On 
the other hand, Aguirre et al. (2004) find little evidence for 
inhomogeneity in the [Si/C] they infer from the observed op- 
tical depth ratio rsnv /tciv , implying that spatial variations 
in the UVB spectral shape are small. There is also some ev- 
idence from metal line analyses for a hardening of the UVB 
spectrum at z ~ 3 (Savaglio et al. 1997; Songaila 1998, 2005; 
Vladilo et al. 2003; Agafonova et al. 2005, 2007), perhaps 
due to the onset of He n reionisation. However, other studies 
are consistent with no evolution at the same redshifts (Kim 
et al. 2002; Boksenberg et al. 2003; Aguirre et al. 2004). 
Quantifying the effect of a spatially inhomogenous UVB on 
intergalactic metal absorption lines is thus desirable, both 
as a potential probe of He n reionisation and for examining 
the systematic bias, if any, spatially uniform ionisation cor- 
rections impart to constraints on the distribution of metals 
in the IGM. 

Analytical arguments have previously highlighted the 
importance of small scale fluctuations in the hydrogen ionis- 
ing radiation field for metals associated with high H i column 
density systems (Miralda-Escude 2005; Schaye 2006). How- 
ever, there has been no detailed examination of the impact 
of a spatially inhomogeneous UVB at E > 4 Ryd on metals 



1 There is some evidence for a reduction in the number of 
C IV absorption systems in close proximity (< 0.3 Mpc) to quasars 
at 1.6 < z < 4 (Wild et al. 2008, see also Fox et al. 2008; Tytler 
et al. 2009). This line-of-sight proximity effect arises because of 
the 1 / R 2 dependence of the ionisation radiation intensity around 
the quasar when the IGM is optically thin, and it will manifest it- 
self even after He II reionisation has long completed. This effect is 
slightly different to the large scale fluctuations in the UVB spec- 
tral shape discussed here, which are expected to be important 
toward the tail-end of He II reionisation when the He II ionising 
photon mean free path is small. 



associated with lower H i column densities. A self-consistent 
treatment would require multi-frequency cosmological radia- 
tive transfer in a large ( > (200 Mpc) 3 ) volume coupled to a 
detailed chemodynamical model. In contrast, current state- 
of-the-art hydrodynamical simulations of IGM metal enrich- 
ment generally assume a spatially uniform UVB (e.g. Tescari 
et al. 2010; Shen et al. 2010; Cen & Chisari 2010; Wiersma 
et al. 2010, but see Oppenheimer et al. 2009) while large cos- 
mological radiative transfer simulations of He n reionisation 
are still restricted to a treatment of hydrogen and helium 
only (Paschos et al. 2007; McQuinn et al. 2009). 

In this paper we instead use a comparatively simple 
model which includes most of the relevant aspects of the spa- 
tially inhomogeneous UVB expected towards the tail-end of 
Hen reionisation (e.g. Fardal et al. 1998; Bolton et al. 2006; 
Furlanetto 2009a). We combine this model with a large hy- 
drodynamical simulation of the IGM to explore the impact 
of spatial fluctuations in the UVB spectral shape on inter- 
galactic carbon and silicon at z = 3. In particular, we focus 
on the ratio of Si iv to C iv absorption, which is sensitive 
to the shape of the UVB either side of the He n ionisation 
edge (Songaila et al. 1995; Savaglio et al. 1997; Giroux & 
Shull 1997). We describe our hydrodynamical simulation and 
model for spatial fluctuations in the UVB spectral shape in 
section 2. Section 3 describes the four toy UVB models we 
use in this work. In section 4 we proceed to examine the 
impact of spatial fluctuations in the UVB spectral shape on 
synthetic absorption spectra constructed from our simula- 
tions. We compare our results to observational data in sec- 
tion 5, and conclude in section 6. All atomic data is taken 
from Morton (2003), solar abundances are from Asplund 
et al. (2009) and all distances are given in comoving units. 



2 NUMERICAL MODEL 

2.1 Hydrodynamical simulation of the IGM 

We use the R4 cosmological hydrodynamical simulation de- 
scribed in Becker et al. (2011) to model the IGM density 
field. The simulation was performed in a 40/i -1 Mpc box 
with 2 x 512 3 particles using the parallel Tree-SPH code 
GADGET-3 (Springel 2005). The cosmological parameters 
are fl m = 0.26, fi A = 0.74, Q, h h 2 = 0.023, h = 0.72, 
as = 0.80, n s = 0.96 (e.g. Komatsu et al. 2009; Reichardt 
et al. 2009). We use a snapshot drawn from the simulation 
at z = 2.976. 



2.2 A simple model for spatial fluctuations in the 
UV background spectral shape 

Our model for spatial fluctuations in the UVB spectral 
shape is similar to the approach described by Bolton et al. 
(2006) (see also Fardal et al. 1998; Furlanetto 2009a). The 
model applies to fluctuations in the UVB spectrum at the 
tail-end of Hen reionisation only, following the overlap of 
He in regions when the He n ionising photon mean free path 
is set by the abundance of Lyman limit systems (Miralda- 
Escude et al. 2000). There is good evidence to suggest that 
Hen reioniation is indeed completing by z ~ 3 (Shull et al. 
2010; Becker et al. 2011), and this assumption should there- 
fore be reasonable. Note, however, that the model will not 
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apply to the patchy ionisation state of the pre-overlap IGM 
during the early stages of Hen reionisation at z > 3, and 
does not include the effect of detailed, small scale radia- 
tive transfer effects on the IGM (Maselli & Ferrara 2005; 
Tittley & Meiksin 2007). Further modelling will thus be re- 
quired to address the impact of spatial variations in the UVB 
spectral shape on intergalactic metals during the heart of 
Hen reionisation. 

Firstly, to create a volume large enough to contain sev- 
eral He ii ionising photon mean free paths, we stack the hy- 
drodynamical simulation volume around the original 40/i _1 
Mpc box to create a cube 600/i _1 Mpc on each side. We 
identify haloes in the volume using a friends-of-friends halo 
finder. Quasar luminosities are assigned to the haloes by 
Monte Carlo sampling the Hopkins et al. (2007) B-band 
quasar luminosity function at z = 3. The luminosities 
are uniquely mapped to the halo masses by rank order- 
ing both quantities and assigning the brightest quasars to 
the most massive haloes. This approach enables us to dis- 
tribute sources on large scales while following the ionisa- 
tion state of gas in detail in the central 40/i -1 Mpc box. 
However, it will only approximately model the clustering 
properties of the sources. Fortunately, the rarity of bright 
quasars means that Poisson fluctuations rather than clus- 
tering is expected to dominate the resulting UVB fluctu- 
ations (Furlanetto 2009a), so this limitation is unlikely to 
significantly affect our results. 

We require the number of quasars in a volume, V, to 
satisfy 

/oc 
<f>(L B )dL B , (1) 

where 4>(Lb) is the Hopkins et al. (2007) quasar luminos- 
ity function. The lower limit of the integral is L m in = 
10 43 ' 5 ergs -1 , corresponding to M B = -20. The B-band lu- 
minosity of each quasar is converted into a luminosity below 
the H i Lyman limit using a broken power law spectrum (e.g. 
Madau et al. 1999) 

!v-° z (2500 < A < 4600 A), 
i/ -0 ' 8 (1050 < A < 2500 A), (2) 
v- a ° (A < 1050 A). 

We consider two cases for the extreme UV (EUV) spectral 
index, a s . In the first case we assume a constant value of 
a s = 1.5, consistent with the mean obtained by Telfer et al. 
(2002) for radio quiet quasars, (a B ) = 1.57 ±0.17. However, 
a wide range of values are measured for the EUV spectral 
index (e.g. Zheng et al. 1997; Telfer et al. 2002; Scott et al. 
2004) . Consequently, we also consider a variable EUV spec- 
tral index. Scatter in a B is included by Monte-Carlo sam- 
pling a Gaussian with mean 1.5 and standard deviation 0.5 
over the range < a s < 3, similar to the dispersion mea- 
sured by Telfer et al. (2002). 

Once the luminosity and spectral energy distribution for 
each quasar is specified, the specific intensity of the UVB, 
J(r, v) [ergs -1 cm -2 Hz -1 sr _1 ], at frequencies between the 
Hi and Hen ionisation edges, v B \ < v < ^Hcii, is given by 

i— 1 

The sum is made over all quasars in the 600ft -1 Mpc volume, 



Table 1. The four UVB models used in this work. From left to 
right, the columns list the model name, the photon mean free 
path at the Hen ionisation edge, the quasar EUV spectral index 
and the logarithm of the spatially averaged softness parameter, 
log<S(r)> = log[(r H i(r)/r H ell(r)>]. 



Model 


^Hell 


a s 


log(S(r)> 


UVB1 


30 Mpc 


1.5 


2.66 


UVB2 


15 Mpc 


1.5 


3.10 


UVB3 


30 Mpc 


Variable 


2.51 


UVB4 


30 Mpc 


1.5 


2.66 



where | ri — r| is the distance of quasar i from r, which is 
always in the central 40/i -1 Mpc box. Eq. (3) assumes that 
the IGM is optically thin to H i ionising photons; this should 
be a reasonable approximation at z = 3. 

At frequencies above the Hen ionisation edge, v > 
I'Heii, the specific intensity is instead 

1 ^ T .( r . v \ |r,-r| / „ \-303-l) 

i — 1 

where AhcII is the mean free path for photons at the 
Hen ionisation edge and /3 is the power-law index de- 
scribing the slope of the Hen column density distribution, 
f(Nneii) oc N^ n . This expression assumes quasars arc 
point sources surrounded by discrete, Poisson distributed 
absorbers (Faucher-Giguere et al. 2009). We adopt /3 = 1.5, 
following the analogous distribution for H i absorbers (Petit- 
jean et al. 1993; Kim et al. 2002). Note, however, the slope 
of the Hen column density distribution is not constrained 
observationally and may deviate from the Hi distribution 
for He ii absorbers associated with high H i column densi- 
ties, JVhi > 10 16 cm -2 (Fardal et al. 1998; Faucher-Giguere 
et al. 2009). Finally, in this work we model the contribution 
to the UVB from quasars only. The non-thermal emission 
from quasars is thought to dominate the UVB at v>VH_ e u, 
but star forming galaxies will make an increasingly signifi- 
cant contribution to the UVB at lower frequencies at z > 3 
(Bianchi et al. 2001; Haehnelt et al. 2001; Bolton et al. 2005; 
Kirkman et al. 2005; Faucher-Giguere et al. 2008). However, 
the soft spectra of these sources are unlikely to significantly 
contribute to spatial fluctuations in the UVB spectral shape 
at E > 4Ryd. 



3 TOY UV BACKGROUND SPECTRA 

The four toy UVB models used in this work are sum- 
marised in Table 1. The models are constructed by eval- 
uating Eqs. (3) and (4) on a 14 3 grid within the central 
40/i -1 Mpc simulation box. The size of each grid cell, ~ 4 
Mpc, is chosen to match the scales of 4 — 10 Mpc on which 
the He n to Hi optical depth ratio is observed to vary at 
z = 2.4 - 2.9 (Shull et al. 2010). This choice also represents 
the best compromise between resolution and efficiency; the 
construction of the photo-ionisation balance look-up tables 
(see section 4.1) becomes very time consuming if the num- 
ber of grid cells much larger. However, note that the rarest, 
high amplitude fluctuations in the UVB spectral shape on 
small scales will not be fully captured in the simulation. In 
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Figure 1. Left: The spatially averaged specific intensity as a function of photon energy, E = hu, for the four spatially fluctuating UVB 
models summarised in Table 1. Models UVB1 (solid curve) and UVB4 (dot-dashed curve) arc identical at E > 4 Ryd. Models UVB1 
and UVB2 (dotted curve) are the same at E < 4 Ryd. Vertical tick marks correspond to the ionisation potentials of various ions. Right: 
The solid curve displays the spatially averaged specific intensity for model UVB1. The range in specific intensity, J(r,u), encompassing 
95 per cent of all fluctuations from the median is displayed by the dotted curves. For comparison, the dashed line shows the UVB 
background model of Haardt & Madau (2001) for emission from galaxies and quasars at z = 3, renormalised to have specific intensity 
J(ym) = 3.5 x 10~ 22 ergs" 1 cm" 2 Hz" 1 sr" 1 . 



each grid cell the specific intensities are evaluated at 25 dif- 
ferent photon energies spanning the range E = 1 — 100 Ryd. 
Lastly, all four models are renormalised by a factor of 0.74 to 
give identical mean specific intensities at the Hi ionisation 
edge (J(r,z/Hi)) = 3.5 x 10~ 22 erg s" 1 cm" 2 Hz" 1 sr" 1 . This 
renormalisation yields H i photo-ionisation rates of Thi ~ 
1.0 x l(T 12 s _1 (cf. Thi ~ 1.3 x 10 _12 s _1 prior to the 
renormalisation), consistent with observational constraints 
derived from the Lya forest opacity at z — 3 (Bolton et al. 
2005; Faucher-Giguere et al. 2008). 

Models UVB1 and UVB2 are constructed using a sin- 
gle EUV spectra index, a B = 1.5. However, the value of 
the mean free path at the Hen ionisation edge, AhcII, is 
rather uncertain. In this work we adopt a fiducial value 
of Ahcii = 30Mpc at z = 3 (e.g. Fardal et al. 1998; 
Miralda-Escude et al. 2000; Bolton et al. 2006), but the 
exact value will depend on the He n to Hi column density 
ratio, r\ — A<hcIi/A?hi, which exhibits significant fluctua- 
tions approaching z = 3 (Zheng et al. 2004; Shull et al. 
2004, 2010). We thus also consider a smaller mean free 
path for model UVB2, Ahcii = 15Mpc, consistent with 
the constraints recently presented by Dixon & Furlanetto 
(2009) at z ~ 3. Model UVB3 instead uses a variable 
EUV spectral index and Aimi = 30 Mpc. The final model, 
UVB4, is identical to model UVB1, aside from the spec- 
tral shape between 3 and 4 Ryd. Resonant absorption due 
to He ii Lyman series lines will also produce significant spa- 
tial fluctuations in the UVB between He n Lya at 3 Ryd 
and the Hen ionisation edge at 4 Ryd (Madau & Haardt 
2009). Model UVB4 mimics these fluctuations by assuming 
the UVB spectrum computed with Eqs. (3) and (4) follows 
a power law J(r, v) = J(r,i^ ya )(^/n, y a) ? for z/L ya < v < 



J4MI, where v-L ya is the frequency of Hen Lya at 3 Ryd and 
£ ~ 81og[J(r, i/Heii)/J(r, H, ya )]. 

The spatially averaged spectra, (J(r, v)), for all four 
UVB models are displayed in the left hand panel of Fig. 1. 
The step at 4 Ryd in models UVB1, UVB2 and UVB3, and 
between 3 — 4 Ryd for model UVB4, is due to the atten- 
uation of the spectrum by the Hen opacity. The spectra 
recover their intrinsic shape at high energies, where the 
photon mean free path is long and the attenuation of the 
radiation by the intervening Hen absorbers is negligible. 
Model UVB3, which assumes a variable EUV spectral in- 
dex, is much harder than the other models at high energies. 
This model is dominated by the quasars with the hardest 
spectra, a s ~ 0, in the simulation volume, even if these ob- 
jects are not the brightest quasars at lower frequencies. The 
logarithm of the spatially averaged softness parameter for 
each model, defined as the ratio of the Hi to Hen photo- 
ionisation rates, S = Thi/ThcH, are listed in Table 1. Ob- 
servational constraints on the softness parameter at z — 3 
are dominated by the large uncertainty in the He n effective 
optical depth. Bolton et al. (2006) find 2.3 < log S < 3.2 
at z — 3 by comparing hydrodynamical simulations of the 
Hi and Hen Lya forest opacity to observational data. All 
four models are designed to be consistent with these con- 
straints. 

The right hand panel of Fig. 1 again displays model 
UVB1, but the dotted curves now show the range en- 
compassing 95 per cent of all fluctuations around the me- 
dian specific intensity. Fluctuations are largest close to the 
He ii ionisation edge where the attenuation of the radiation 
is strongest, while towards higher frequencies the increasing 
mean free path produces much smaller departures from the 
median. For comparison, the dashed line shows the widely 
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Figure 2. Left: The probability distribution for fluctuations in the specific intensity, J(r, z-'Hcll): relative to the spatially averaged value in 
the simulation volume, ( J(r, fH e ii)), for the four UVB models summarised in Table 1. The distributions are calculated using the specific 
intensity at the Hen ionisation edge, where the fluctuations are largest. Note that the distributions for models UVB1 and UVB4 arc 
identical at this frequency. Right: Cumulative probability distributions for fluctuations in the specific intensity of UVB2 at 4 Ryd. The four 
curves show the cumulative distributions obtained in pixel subsets corresponding to overdensities at log A = — 1, 0, 1 and 2. 



used UVB model of Haardt & Madau (2001) (HM01) for 
emission from galaxies and quasars at z = 3. The spec- 
trum has been renormalised to match the specific intensity 
of UVB1 at 1 Ryd. The softness parameter of the HM01 
spectrum is logS 1 = 2.42, around 0.2 dex lower than UVB1. 
Although the features of the toy UVB spectra presented here 
are broadly similar to the more detailed model of HM01, 
there are some important differences. In particular, our toy 
models do not include more complex effects such as recombi- 
nation emission (Haardt & Madau 1996; Fardal et al. 1998; 
Faucher-Giguere et al. 2009) or sawtooth absorption by the 
Hen Lyman series (Madau & Haardt 2009) which are im- 
portant to consider when performing detailed modelling of 
observational data. 

The fluctuations in the specific intensity are quan- 
tified in more detail in Fig. 2 for all four UVB mod- 
els. The left panel displays the probability distribution 
for the specific intensity relative to its spatially averaged 
value, log[J(r, i/Heii)/(J(r, ^Heii))], at the Hen ionisation 
edge. The distribution becomes broader for a smaller AhcII 
(UVB2) as the number of sources within one mean free path 
is reduced. Somewhat counter-intuitively, the variable EUV 
spectral index model (UVB3) produces fewer of the large, 
rare fluctuations from the mean than UVB1, which assumes 
a constant a s but is otherwise identical. This is because the 
three quasars which lie closest to the centre of our simula- 
tion volume (and hence dominate the rare, high amplitude 
fluctuations) have randomly assigned EUV spectral indices 
which conspire to reduce fluctuations in the specific inten- 
sity at 4 Ryd. Recall the quasar luminosities are normalised 
in the B-band, and in this instance the EUV spectra for the 
three closest quasars all have a s > 1.5. This produces lower 
specific intensities at the He n ionisation edge compared to 
model UVB1 with a B = 1.5, and hence fewer large fluctua- 
tions from the mean. 



Lastly, the right panel of Fig. 2 displays cumulative 
distributions for log[J(r, vn e n)/{J{r, ^Heii))] obtained from 
model UVB2. The four separate curves show the cumulative 
distributions for subsets of pixels corresponding to regions in 
the IGM density field with overdensities at log A = — 1, 0, 1 
and 2, respectively. The lack of a clear correlation with den- 
sity 2 is mainly due the low number density of bright quasars, 
resulting in a radiation field which varies on scales much 
larger than the Jeans scale (see the discussions in e.g. Mc- 
Quinn et al. 2009; Furlanetto 2009b). Note, however, that 
fluctuations on scales < 4Mpc (typically those associated 
with the radiation field in close proximity to the quasars) will 
be smoothed over, but these fluctuations are also rather rare. 
On the other hand, if additional physics such as small scale 
radiative transfer effects are important, or if Hen ionising 
sources are far more numerous (e.g. galaxies), the ionising 
radiation field may be more closely correlated with the IGM 
density on small scales. 



2 Bolton et al. (2006) found that the He II to Hi column den- 
sity ratio, 77 = Nn e ii/Nni, obtained from a Voigt profile analysis 
of synthetic absorption spectra was systematically higher in re- 
gions where tjji > 0.05 compared to regions with lower Hi optical 
depths (and hence densities). They suggested this implied the 
fluctuating UVB model they used was harder in higher density 
regions. However, this interpretation is in disagreement with our 
findings here. This earlier result is likely spurious, and probably 
results from the assumptions of pure turbulent broadening in the 
Voigt profile analysis (Fechner & Reimers 2007) and a spatially 
uniform H I photo-ionisation rate. 
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Figure 3. Upper left: The cumulative probability distribution for thcIi/thi for all pixels in the sight-lines drawn from our hydrodynamical 
simulation. The curves correspond to models UVB1 (solid), UVB2 (dotted), UVB3 (dashed) and UVB4 (dot-dashed). Thick curves are 
obtained from spectra constructed using the spatially fluctuating UVB. Thin curves correspond to the optical depth ratios obtained using 
the spatially averaged UVB spectrum. Upper right: As for the upper left panel, but now for pixels with 0.02 < THi.Hell < 5 on ly- Models 
UVB1 and UVB4 are indistinguishable in both of these panels. Lower left: The cumulative probability distribution for TgnvA'Crv f° r a U 
pixels. Lower right: As in the lower left panel, but now for pixels with 0.02 < Tciv.Silv < 5 only. The thin and thick curves are almost 
indistinguishable in the lower panels. Note also the scales on the horizontal axes of the lower left and right hand panels arc significantly 
different. 



4 THE IMPACT OF UVB FLUCTUATIONS ON 
SIMULATED ABSORPTION SPECTRA 

4.1 Construction of synthetic absorption spectra 

We now turn to examine the impact of our spatially in- 
homogeneous UVB models on synthetic quasar absorption 
spectra. We use a simple a posteriori approach to includ- 
ing metals within our hydrodynamical simulation. We do 
not model the production and dispersal of metals in de- 
tail (see e.g. Theuns et al. 2002a; Cen et al. 2005; Op- 
penheimer & Dave 2006; Shen et al. 2010; Wiersma et al. 
2010; Tescari et al. 2010). We instead follow carbon and 
silicon abundances only, and assume the metallicity of the 
IGM traces the gas density distribution. Schaye et al. 
(2003) measure a best fit median metallicity of [C/H] = 
-3.47tooe + 0.65j;o 12 (log A - °- 5 )> with a lognormal scat- 
ter of <t([C/H]) = 0.76lg;g|-0.23tg;g?(logA-0.5) at z = 3. 
Following Schaye et al. (2003), we introduce the lognormal 
scatter to the median as [C/H] = -3.47+0.65(log A-0.5)+s 
by dividing the 40/i -1 Mpc simulation volume into 32 3 cubes 
and adding a different scatter, s, in each sub-volume. The 
scatter, s, is obtained by randomly sampling a normal dis- 
tribution in [C/H] with zero mean and standard deviation 
a = 0.76. We further assume the silicon abundance rela- 
tive to carbon is [Si/C] = 0.77 (Aguirre et al. 2004). Note, 



however, these metallicity constraints are derived from a sta- 
tistical analysis of C iv and Si iv absorption lines which use 
an ionisation correction based on the HM01 UVB model. 
Different ionisation corrections obtained using significantly 
softer or harder UVB spectra will alter these constraints 
considerably. 

The equilibrium ionisation balance of the metals in our 
hydrodynamical simulation is calculated using the photo- 
ionisation code Cloudy (version 08.00), last described by 
Ferland et al. (1998). For each of our four spatially inho- 
mogeneous UVB models, we create a five dimensional look- 
up table listing the ionisation fractions of hydrogen, helium, 
carbon and silicon as a function of gas temperature, T, hy- 
drogen number density, tih, and position on the 14 3 UVB 
grid within our AQhT 1 Mpc simulation box. Synthetic spec- 
tra are obtained from the hydrodynamical simulations by 
randomly selecting 1000 sight-lines from the simulation and 
performing an interpolation on the particle data weighted 
by the smoothing kernel {e.g. Theuns et al. 1998). Each line 
of sight is drawn parallel to the box boundaries and has 2048 
pixels, each of which has an associated gas density, tempera- 
ture, peculiar velocity and ionisation fraction. The ionisation 
fractions for H i, He n, C iv and Si iv are obtained by linearly 
interpolating the five dimensional look-up table. Absorption 
spectra are constructed for Hi Lya (A1216), Ly/3 (A1026), 
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Ly-y (A973), LyS (A950), Hen Lya (A304), Civ (AA1548, 
1551) and Siiv (AA1394, 1403). 

4.2 Optical depth ratios 

The ratio of the optical depths for two different species is 
sensitive to fluctuations in the spectral shape of the ionis- 
ing background between their respective ionisation poten- 
tials. A ratio also has the advantage of being less suscepti- 
ble to fluctuations in the IGM density field (Worseck et al. 
2007; Furlanetto & Lidz 2010). We examine timi/thi for the 
Lya transitions and tshv/tciv for A1548 and A1394 in our 
synthetic spectra. Both of these ratios have been widely used 
as probes of Hen reionisation and the UVB spectral shape 
at 2 ~ 3 (e.g. Songaila 1998; Schaye et al. 2003; Worseck 
et al. 2007; Shull et al. 2010). Note the latter will also be 
sensitive to fluctuations in the relative abundance of carbon 
and silicon. 



4-2.1 The thoIi/thi ratio 

It is instructive to first revisit the effect of UVB fluctuations 
on the He n to H i optical depth ratio (but see also Fardal 
et al. 1998; Maselli & Ferrara 2005; Bolton et al. 2006; Titt- 
ley & Meiksin 2007; McQuinn et al. 2009; Furlanetto 2009a). 
The upper left panel of Fig. 3 shows the cumulative distri- 
bution for timi/thi for all pixels in our synthetic spectra. 
This includes pixel optical depths that are in practice un- 
observable, either because of saturated absorption or noise. 
The thick curves correspond to the optical depth ratios for 
the four spatially fluctuating UVB models listed in Table 1. 
Models UVB1 and UVB4 are almost identical; neither Hi or 
He ii arc significantly ionised by photons with E = 3 — 4Ryd. 
Models with progressively harder spectra produce lower me- 
dian values for the He n to Hi ratio. The He n fraction de- 
creases as the specific intensity at E > 4Ryd increases; since 
the UVB is almost uniform at 1 Ryd in all our models (see 
the right panel in Fig. 1), fluctuations in the specific in- 
tensity at the Hi ionisation edge will have little impact on 
thcIi/thi- Note also that the smaller mean free path used 
for UVB2 produces a much wider range of thcIi/thi val- 
ues due to the wider range of fluctuations produced at the 
He ii ionisation edge. 

We may assess the extent to which thcIi/thi is influ- 
enced by fluctuations in the UVB spectral shape by com- 
parison to the ratio obtained for a spatially uniform UVB. 
The thin curves display the optical depth ratios obtained 
using the spatially averaged (i.e. spatially uniform) spec- 
trum for each UVB model. The distributions for models 
including UVB fluctuations are indeed broader, especially 
toward higher values of thcIi/thi- However, all four distri- 
butions still extend to low thcIi/thi in the spatially uniform 
case. This behaviour is almost entirely due to variations 
in thcIi/thi introduced by thermal broadening; the ther- 

— 1/2 

mal width of an absorption line for species i is 6i oc m ; , 
where mi is the mass of the particle. Fluctuations in the op- 
tical depth ratio therefore also originate from the wings of 
absorption lines which are predominantly thermally rather 
than turbulently broadened (Fechner & Reimers 2007) . 

The upper right panel of Fig. 3 instead shows only the 
subset of pixels which have 0.02 < thi.imi < 5, roughly 



corresponding to range over which the optical depths may 
be reliably measured in observational data (e.g. Shull et al. 
2010). These pixels correspond to 27.4, 4.1, 39.1 and 27.1 
per cent of the total for models UVB1, UVB2, UVB3 and 
UVB4, respectively. The median t Heii /t hi for UVB2 is now 
significantly smaller. This is because the small number of 
optical depths with 0.02 < thcII < 5 selected for UVB2 tend 
to probe regions where the Hen fraction is lowest. 

Recent observations of the Hen Lya opacity in the 
spectrum of the z ~ 2.9 quasar HE2345— 4342, made us- 
ing the Cosmic Origins Spectrograph on the Hubble Space 
Telescope, reliably measure fluctuations in the optical depth 
ratio spanning the range thoIi/thi ~ 1.25—125 at 2.75 < z < 
2.89 (Shull et al. 2010). In comparison, model UVB2 exhibits 
thch/thi ^ 0.8 - 248 (15 - 188 for 95 per cent around the 
median), while UVB1 and UVB4 have thoIi/thi ^ 0.4 - 77 
(16 - 65) and UVB3 has t H cIi/thi ^ 0.3 - 58 (13 - 46). 
We do not perform a more detailed comparison here, but as 
previously noted these results strongly suggest the observed 
t Heii/ thi fluctuations are largely attributable spatial vari- 
ations in the He n ionising background. These fluctuations 
arise primarily from the small number of quasars within the 
short Hen ionising photon mean free path, AhoII ~ 15 — 30 
Mpc, expected towards the tail-end of He n reionisation. 

4-2.2 The tsuv/tciv ratio 

The results for the tsuv/tciv ratio are shown in the lower 
panels of Fig. 3. The ionisation edges for Si iv and C iv are 
at 3.32 and 4.74 Ryd respectively; the ratio of these optical 
depths is thus sensitive to the shape of the UVB either side 
of the Hen ionisation edge at 4 Ryd (Songaila et al. 1995; 
Giroux & Shull 1997). The left hand panel displays the dis- 
tributions for all pixels in the synthetic spectra. Models with 
progressively harder spectra produce smaller median values 
for the ratio - the distribution for UVB3 in particular ex- 
tends to significantly lower rsiiv/TCiv values. A harder UVB 
reduces both the Si iv and C iv ionisation fractions, but the 
Si iv fraction is lowered more rapidly as the break at 4 Ryd is 
reduced, leading to the smaller values for the rsnv /tciv ratio 
in UVB3. However, in contrast to the case for THeii/THi, the 
distributions for the fluctuating and uniform UVB models 
are almost identical. Only the UVB3 distribution displays 
a small difference for the lowest values of rsnv /tciv- This 
scatter in the tsuv/tciv ratio is instead dominated by the 
different dependences of the C iv and Si iv fractions on gas 
density, and to a much lesser extent variations in the ab- 
sorption profiles in redshift space. 

The lower right hand panel in Fig. 3 displays the subset 
of pixels with 0.02 < rsnv, civ < 5. These pixels represent a 
very small fraction of the total: 0.2, 0.3, 0.1 and 0.2 per cent 
for models UVB1, UVB2, UVB3 and UVB4, respectively, re- 
flecting the much smaller volume filling factor of detectable 
C iv and Si iv absorption in the simulations. The medians 
of all four distributions are now significantly larger; most of 
the discarded pixels have C iv and Si iv optical depths which 
are too low to be observed. The distributions for models 
UVB1, UVB2 and UVB3 are now similar, suggesting that 
the rsiiv/Tciv ratio is a less sensitive probe of the UVB 
spectral shape for these pixels. Finally, all four distributions 
are again virtually indistinguishable from the distributions 
obtained with a spatially uniform UVB. 
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Figure 4. Contour plots of the optical depth weighted temper- 
ature, T T , against the optical depth weighted overdensity, A T , 
for pixels with 0.02 < r < 5 in the synthetic absorption spectra 
constructed using model UVB1. Top left: The temperature den- 
sity plane for H I absorption. Top right: He II absorption. Bottom 
left: C IV absorption. Bottom right: Si IV absorption. The number 
density of pixels increases by 1.0 dex within each contour level. 



4.3 Gas temperature and density 

We may gain insight into this behaviour by first examin- 
ing the typical gas densities and temperatures which pro- 
duce the absorption in our synthetic spectra. Contour plots 
of the optical depth weighted temperature against optical 
depth weighted overdensity are shown in Fig. 4 for spec- 
tra constructed using model UVB1. Clockwise from the up- 
per left, each panel shows the temperature-density plane for 
Hi, Hen, Siiv and Civ optical depths with 0.02 < r < 
5. The number of pixels increases by an order of magni- 
tude within successive contours. It is evident that Hi and 
Hen absorption primarily probes photo-ionised gas around 
mean density and below, with the Hen absorption slightly 
more sensitive to gas in voids (Croft et al. 1997; McQuinn 
2009). In contrast, Civ and especially Siiv tend to probe 
overdense regions in the IGM (e.g. Rauch et al. 1997). Nearly 
all the carbon and silicon absorption is furthermore associ- 
ated with gas at T < 10 5 K. 

The exact temperature and density of regions with mea- 
surable optical depths depends on the UVB model used. 
For example, for UVB3 (the hardest spectrum we con- 
sider) the median thi weighted temperature and overden- 
sity for pixels with 0.02 < thi,h b ii < 5 are T T = 7100K 
and log A T — —0.6. For 0.02 < rciv.suv < 5 the me- 
dian tciv weighted values are instead T T — 34 500 K and 
log A T = 2.0. For UVB2 (the softest spectrum) the cor- 
responding values are TV — 5 600K, logA T = —0.8 for 
0.02 < THi.Heii < 5 and TV = 43 400 K, log A T = 1.6 for 
0.02 < tciv, Siiv < 5. Nevertheless, it is clear the thi/timi 
and tsuv/tciv ratios probe rather different temperature and 
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Figure 5. The ratio of the Si IV to C IV ionisation fractions, 
/siiv//civ> f° r * ne four UVB models considered in this work. 
The ratios are plotted as a function of J (ism) / J (vn c u) , the ratio 
of specific intensities at the H I and He II ionisation edges (see 
text for details). In each panel, /siiv//civ is shown at four 
different gas densities, njj = 10 — 2 cm - 3 (dot-dashed curves), 
njj = 10 — 3 cm -3 (solid curves), njj = 10~ 4 cm~ 3 (dotted curves) 
and nn = 10 — 5 cm - 3 (dashed curves). The thick curves assume 
a gas temperature of T = 32 500 K while the thin curves cor- 
respond to T = 65 000 K. The horizontal extent of the shaded 
regions display the range encompassing 95 per cent of all fluctu- 
ations around the median J (vm) / J (vn a ii) f° r the spatially inho- 
mogeneous UVB models. Upper left: UVB1. Upper right: UVB2. 
Lower left: UVB3. Lower right: UVB4. 



density regimes in the IGM. Furthermore, since the majority 
of the C iv and Si iv absorption in our models corresponds 
to gas with T < 10 5 K, collisional ionisation is unlikely to ex- 
plain the similarity between the rsiiv/Tciv distribution for 
the spatially uniform and fluctuating cases shown in Fig. 3. 
Note, however, additional heating from feedback is not in- 
cluded in our hydrodynamical simulation, and we thus likely 
underpredict the amount of C iv and Si iv absorption due 
to collisionally ionised gas. However, because fluctuations 
in the photo-ionising background will be less important for 
the rsiiv/Tciv ratio if a larger fraction of these elements 
are in a hot, predominantly collisionally ionised phase with 
log A > 2, this is unlikely to significantly alter our conclu- 
sions regarding the importance of inhomogeneities in the 
UVB spectral shape. 



4.4 The UVB spectral shape 

With the typical gas densities and temperatures of the ab- 
sorption in hand, we may now examine the effect of varying 
the spectral shape of the UVB in more detail. It is instructive 
to simply rescale the spatially averaged UVB spectrum for 
each model by a constant factor above the Hen ionisation 
edge: 
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Figure 6. Individual ionisation fractions for several different ion- 
isation states of carbon (upper panels) and silicon (lower panels 
for model UVB2 (left column) and UVB4 (right column). The 
ionisation fractions are plotted as a function of J (vm) / J (vn e ii) 
for gas nn = 10~ 3 cm~ 3 and T = 32 500 K. The horizontal ex- 
tent of the shaded regions display the range encompassing 95 per 
cent of all fluctuations around the median J (vm) / J X^Hcll) f° r 
the spatially inhomogencous UVB models. 
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(5) 



where \ is a dimensionless constant. For model UVB4, this 
rescaling also modifies the spectral shape between E — 
3 — 4 Ryd. Smaller (larger) values of x thus produce softer 
(harder) UVB spectra and a stronger (weaker) break in the 
spectrum at 4 Ryd. 

Fig. 5 displays the Siiv to Civ ionisation fraction, 
/siiv//civ, for the four UVB models considered in this work. 
Ignoring redshift space distortions, this is related to the op- 
tical depth ratio by 



TSilV Q~SiIV ?tSi /silV 

tciv crciv nc /civ 



1 j f SilV 10 [Si/C]-0.77 

/civ 



(6) 



where o\ and rn are the absorption cross-section and num- 
ber density for each species. The ratios are plotted as a 
function of the ratio of specific intensities at the Hi and 
He ii ionisation edges, J (um) / J (unai) , following the above 
rescaling. The thick curves are computed using Cloudy for 
gas with T — 32 500 K, similar to the temperature of gas 
responsible for C iv and Si iv absorption in our simulated 
spectra. In each panel, the Siiv to Civ ratio is shown at 
four different gas densities, n H = 10 -2 , 10~ 3 , 10~ 4 and 
10 -5 cm -3 . These have overdensities at z = 3 of log A = 
2.91, 1.91, 0.91 and —0.09, respectively. For comparison, the 
median C iv optical depth weighted overdensity probed by 
Siiv and C iv absorption with 0.02 < rciv.Siiv < 5 for 
UVB1 is log A T = 1.7. 



The horizontal extent of the shaded regions in Fig. 5 
correspond to the range encompassed by 95 per cent of all 
fluctuations from the median J (vni) / J (vncii) in the spa- 
tially fluctuating UVB models. This provides a rough guide 
to the range in fsnv/fciv present in the synthetic spec- 
tra as a function of density at fixed temperature. However, 
there are three important points to keep in mind. Firstly, 
variations in the gas temperature 3 will produce additional 
fluctuations in fsnv/fciv at fixed density. This is illus- 
trated by the thin curves in Fig. 5, which are computed 
for gas with T — 65 000 K; this temperature lies toward 
the upper range of the scatter observed in Fig. 4. Secondly, 
only values of —1.6 < log(Tsiiv/fciv) < 1 (lower right panel, 
Fig. 3) are readily observable in our synthetic spectra, which 
from Eq. (6) corresponds to — 1.8 < log(/siiv//civ) 0.8. 
Thirdly, the range of fsnv/fciv fluctuations will in practice 
be smaller than indicated by the range of J(t'Hi)/ J(^Hcii) 
shown by the shading in Fig. 5. Fluctuations in the UVB 
spectral shape are maximised at 4 Ryd where the mean free 
path is shortest (e.g. Fig. 1), whereas Civ and Siiv have 
ionisation potentials above and below 4 Ryd, respectively. 
Furthermore, the fluctuations do not have an equal proba- 
bility of occurring in this range, and have a higher proba- 
bility of lying close to the median (e.g. Fig. 2). The shaded 
regions thus correspond to an upper limit in the expected 
/siiv//civ variation for each model. 

With these points in mind, the /snv//civ ratio never- 
theless shows only modest variation (< 1.0 dex) over the 
shaded range in Fig. 5 for nn = 10~ 3 cm~ 3 , with the largest 
variations occurring for UVB2 and UVB4. Towards larger 
J (vm) / J (vKeii) (i.e. softer spectra), fsnv/fciv noticeably 
flattens for all models, and becomes almost constant for 
UVB1, UVB2 and UVB3. Larger variations in fsnv/fciv do 
occur at lower densities, especially for model UVB3 where 
the Siiv to Civ ratio changes by ~ 2.2 dex within the 
shaded range for «h = 10~ 5 cm" 3 . The Siiv to Civ optical 
depth ratio for all pixels in Fig. 3 is slightly broader than the 
uniform case for UVB3 at low values as a consequence. How- 
ever, from Eq. (6) we infer that C iv and Si iv absorption 
from such low density gas is not detectable in the synthetic 
spectra. 

In Fig. 6 the individual carbon (upper panels) and sil- 
icon (lower panels) ionisation fractions for several differ- 
ent ions are shown against J(i'hi)/</(j4mi). The curves are 
again computed assuming densities and temperatures typi- 
cal of the C iv and Si iv absorption in our simulated spectra, 
nn = 10~ 3 cm" 3 and T = 32 500 K. The left panels display 
the results for UVB2, while UVB4 is shown on the right. 
Models UVB1 and UVB3 display similar behaviour to UVB2 
within the equivalent ranges for J(^Hi)/J(^Heii). Note the 
C iv and Si iv fractions for model UVB2 are largely insensi- 
tive to the UVB spectral shape for log[J(z/Hi)/J(z/H e ii)] > 2 
and 3, respectively. This is because the UVB is very soft, 



3 We do not include the effect of large scale ( > 50 Mpc) spatial 
fluctuations in the temperature of the low density IGM expected 
during He II reionisation (Theuns ct al. 2002b; McQuinn et al. 

2009) . However, these fluctuations are expected to have a signifi- 
cant impact on large scale, transverse correlations in the Lya for- 
est only; their signatures are much more difficult to detect along 
individual lincs-of-sight and on smaller scales (McQuinn et al. 

2010) . 
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resulting in photo-ionisation rates that are too low to sig- 
nificantly change the C iv and Si iv fractions in the higher 
density gas responsible for most of the absorption. On the 
other hand, the Civ fraction for model UVB4 is more sen- 
sitive to an increase in the strength of the 3 — 4 Ryd break. 
The break at 3 Ryd in this model produces a larger pro- 
portion of carbon in the form of C in relative to C iv for 
increasing J (vm) / J (vhcIi) (Agafonova et al. 2007; Madau 
& Haardt 2009; Vasiliev et al. 2010). However, Siiv again 
flattens towards larger values for J (vhi) / J (vkcu) , moderat- 
ing the change in the Si iv to C iv ratio. Furthermore, as the 
C iv fraction drops rapidly at \og[J(vm) / J(vndi)] > 2 it will 
become progressively more difficult to detect the weakening 
C iv absorption. 

Consequently, as the Hen ionising photon mean free 
path becomes smaller toward higher redshift, the average 
UVB spectral shape softens and variations in the observed 
Si iv to C iv ratio do not become significantly more pro- 
nounced as a result of spatial fluctuations in the UVB spec- 
tral shape. The Civ and Siiv absorption systems typically 
originate from overdense regions with log A ~ 2 in our simu- 
lations, where the Si iv and C iv fractions change less rapidly 
with increasing J(vhi)/J(vh c ii) compared to lower density. 
It is this behaviour, combined with the smaller fluctuations 
expected in the UVB spectral shape at frequencies above 
and below the Hen ionisation edge, which explain the sim- 
ilarity between the Si iv to C iv optical depth ratios for the 
uniform and fluctuating models in Fig. 3. This is in strik- 
ing contrast to the He n to Hi ratio, where the same spa- 
tial variations in the UVB spectral shape significantly in- 
crease fluctuations in thcIi/thi from lower density gas as 
the mean free path is lowered (Fardal et al. 1998; Bolton 
et al. 2006; Furlanetto 2009a). The absence of any observa- 
tional evidence for fluctuations in the Tsiiv/i"civ ratio (or 
the inferred silicon to carbon ratio) is thus not necessarily 
indicative of a spatially uniform UVB at z ~ 3. On the other 
hand, slightly larger tsuv/tciv fluctuations are expected at 
lower densities, although the weaker absorption from these 
regions is much more difficult to detect. We now examine 
whether UVB fluctuations can be detected statistically in 
the observational data, and whether or not they contribute 
to the observed scatter in the IGM metallicity (Rauch et al. 
1997; Schaye et al. 2003; Simcoe et al. 2004). 



5 PIXEL OPTICAL DEPTH ANALYSIS 
5.1 Method 

The pixel optical depth (POD) procedure provides a pow- 
erful tool for statistically analysing metal absorption at low 
gas densities (e.g. Songaila 1998; Ellison et al. 2000; Pieri & 
Haehnelt 2004). Briefly, the procedure involves recovering 
optical depths for a base transition (e.g. Hi Lya absorp- 
tion) on a pixel- by- pixel basis and pairing these with metal 
optical depths recovered at the same redshift. Given suit- 
able calibration, typically achieved using a cosmological hy- 
drodynamical simulation, information may then be derived 
on the abundance and distribution of metals in the IGM. 
The C iv POD analysis performed by Schaye et al. (2003) 
found evidence for scatter in [C/H] at fixed thi, described 
by a normal distribution in [C/H] with standard deviation 



tr([C/H]) = 0.76 - 0.23(log A — 0.5) at z = 3. Scannapieco 
et al. (2006) and Pieri et al. (2006) noted this scatter may 
be attributable to spatial variations in the metallicity of the 
IGM, although they did not rule out the possibility that 
spatial fluctuations in the UVB spectral shape (and hence 
the assumption of a uniform ionisation correction) may also 
play a role. On the other hand, Aguirre et al. (2004) per- 
formed a POD analysis for Siiv and found no evidence for 
additional scatter in the tsuv/tciv ratio beyond the afore- 
mentioned scatter in the metallicity. The study concluded 
that a uniform [Si/C] is favoured and that inhomogeneities 
in the spectral shape of the UVB are small. However, none 
of these studies modelled the effect of UVB fluctuations on 
the ionisation balance for carbon and silicon in detail. 

We perform a POD analysis on our synthetic spectra 
using the procedure outlined by Aguirre et al. (2002). We 
consider the optical depths for Hi Lya through to LyS as 
well as C iv and Si iv absorption. For the POD analysis we 
construct absorption spectra which resemble high resolution, 
high signal-to-noise observational data (e.g. Schaye et al. 
2003). The synthetic spectra are convolved with a Gaussian 
instrument profile of width 7kms~ 1 , resampled onto pix- 
els of width 3.1 kms -1 and Gaussian distributed noise with 
S/N — 100 is added. Continuum fitting errors are mimicked 
by performing an iterative continuum correction to the syn- 
thetic spectra. We compute the median transmitted flux in 
each synthetic line-of-sight and deselect all pixels below la 
of this value, where a is the rms noise amplitude in each 
pixel. The median flux is computed again for the remaining 
pixels, and the procedure is then continued until convergence 
is reached. 

5.2 Comparison to observational data 

The results of our POD analysis are presented in Fig. 7. 
The data points with la error bars in the left panel cor- 
respond to the Schaye et al. (2003) observational measure- 
ments for Q1422+230. These data have a median absorption 
redshift of z = 3.225, at slightly higher redshift than our 
synthetic absorption spectra. From top to bottom, the three 
sets of data points correspond to the 84th, 69th and 50th 
percentiles of the recovered C iv optical depths; the 84th and 
69th percentiles have been offset by +1.5 dex and +0.75 dex 
for clarity. The curves display the recovered C iv pixel opti- 
cal depth against H i optical depth for spectra constructed 
using model UVB4, which we find has the largest impact 
on the recovered optical depths. The solid curves are ob- 
tained from the spatially inhomogeneous UVB model, while 
the dotted curves corresponds to the optical depths recov- 
ered from the spatially averaged spectrum. The curves are 
almost indistinguishable, indicating that the predicted spa- 
tial inhomogeneities in the UVB spectral shape approaching 
He ii reionisation will not significantly impact on the ionisa- 
tion correction for C iv. 

The dashed and dot-dashed curves in Fig. 7 use the spa- 
tially uniform UVB model, but now also include different 
assumptions for the lognormal scatter in the IGM metallic- 
ity at fixed density. The dashed curves assume zero scatter, 
a — 0, while the dot-dashed curves correspond to a = 1.52. 
The impact on the recovered optical depths is especially 
prominent for the higher percentiles, suggesting that spa- 
tial variations in metallicity rather than the UVB spectral 
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Figure 7. Left panel: Comparison of the C IV pixel optical depths as a function of H I optical depth, obtained from model UVB4 to the 
Schaye et al. (2003) observational measurements for Q1422+230. From top to bottom, the data points correspond to the 84th, 69th and 
50th percentiles of the recovered C IV optical depths. The 84th and 69th percentiles have been offset by +1.5 dex and +0.75 dex for clarity. 
The solid curves correspond to the pixel optical depths recovered using the spatially inhomogeneous UVB model, while the dotted curves 
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drawn from a distribution with <r = 1.52, twice the fiducial value of a = 0.76 . Right panel: The recovered Si IV optical depth as a function 
of Civ optical depth for UVB4. The simulations arc compared to the observational data of Aguirre et al. (2004) for Q1422+230. The 84th 
and 69th percentiles are again offset by +1.5 dex and +0.75 dex for clarity, and the different curves are as described for the left hand panel. 



shape will dominate any scatter. As noted by Schaye et al. 
(2003), the fiducial model with a = 0.76 is in somewhat 
better agreement with the higher percentiles. 

In the right hand panel we compare the same model 
to the POD measurements for Si iv and C iv presented by 
Aguirre et al. (2004) for the same quasar. Again, the re- 
covered optical depths are very similar for the spatially in- 
homogeneous and uniform UVB models, and scatter in the 
metallicity has a much larger impact on the recovered op- 
tical depths. The fiducial model with a — 0.76 is again in 
better agreement with the higher percentiles, although the 
agreement is poorer for the median. Note also, in contrast 
to the Civ and Hi pixel optical depths, lognormal scatter 
in the metallicity at fixed density lowers each percentile for 
rsiiv as a function of tciv- The explanation for this be- 
haviour is that scatter in metallicity is added to both of the 
pixel optical depths, in contrast to the left hand panel of 
Fig. 7. Adding this scatter enables the detection of both 
C iv and Si iv absorption in more pixels at lower densities 
due to their higher metallicities; these pixels are otherwise 
hidden in the flat (noise and continuum error dominated) 
part of the correlation at logrciv < —1.5. However, for 
this UVB model the Si iv fraction (and hence optical depth) 
decreases more rapidly with decreasing gas density rela- 
tive to the C iv fraction. Consequently, a greater fraction of 
C iv optical depths with low rsiiv are now detectable, leading 
to the lowering of the rsnv percentiles observed in Fig. 7. 

Finally, we caution the reader not to take the differ- 
ences between the observational data and simulations too 
seriously; the Schaye et al. (2003) and Aguirre et al. (2004) 



metallicities we have used to construct our spectra are de- 
rived under the assumption of a different UVB model. There 
is therefore no reason to expect our model to match the data 
exactly using these metallicities; the comparison here is in- 
stead largely illustrative and one may always fine tune the 
metallicity to enable a given UVB model to match the data. 
The key point is that we clearly expect variations in the 
spatial distribution of metals to dominate any scatter in the 
C iv and Si iv optical depths at fixed density, even in the 
presence of the spatially inhomogeneous UVB expected at 
the tail-end of Hen reionisation. 



6 CONCLUSIONS 

We use a large hydrodynamical simulation of the IGM com- 
bined with a toy model for spatial inhomogeneities in the 
UVB spectral shape to investigate the impact of spectral 
hardness fluctuations on the ionisation balance of intergalac- 
tic carbon and silicon at z ~ 3. We construct synthetic 
quasar absorption spectra from the simulations, assuming 
the metallicity of the IGM traces the underlying gas density 
(Schaye et al. 2003; Aguirre et al. 2004). We carefully ex- 
amine the impact of the UVB fluctuations on the Siiv and 
C iv optical depths. Four different spatially inhomogeneous 
UVB models which employ a variety of different assumptions 
for the UVB spectral shape are considered in our analysis. 
We reconfirm that fluctuations in the UVB spectral shape 
expected at the tail-end of He n reionisation have a signifi- 
cant impact on the He n to H i optical depth ratio (see also 
Fardal et al. 1998; Bolton et al. 2006; Furlanetto 2009a). 
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However, some of the lowest values for this ratio also re- 
sult from lines which are predominantly thermally broad- 
ened (Fechner & Reimers 2007). On the other hand, while 
the Si iv to C iv ratio is indeed sensitive to the average spec- 
tral shape of the UVB, we find the predicted fluctuations 
have little impact on tsuv/tciv measured in our synthetic 
spectra. 

The majority of the detectable Siiv and C iv absorption 
in our synthetic spectra originates from regions with log A ~ 
1.5 — 2 and T ~ 35 000 K. These absorbers are predominantly 
photo-ionised in our simulations, and we may thus exclude 
collisional ionisation as a possible explanation for the small 
impact of UVB fluctuations on our simulated spectra. In- 
stead, we find the ratio of observable Si iv to C iv optical 
depths varies relatively little considering the wide range of 
fluctuations in the UVB spectral shape. This is in part be- 
cause of the longer mean free path for photons above and 
below the He n ionisation edge, which results in smaller fluc- 
tuations in the UVB spectral shape at these frequencies. 
However, as fluctuations in the UVB spectral shape become 
larger as the Hen opacity increases, the spatially averaged 
UVB spectral shape becomes softer; UVB models which pro- 
duce the largest fluctuations in thcIi/thi also have photo- 
ionisation rates which are too low to have a significant im- 
pact on the observed Si iv to C iv ratio. At lower gas den- 
sities, or for UVB models which predict a larger fraction 
of Cm relative to Civ (e.g. for Hen Lyman series absorp- 
tion), the expected variation in the Siiv to Civ ratio can 
be slightly larger, but it is more difficult to detect due to 
the correspondingly smaller gas densities and/or fraction of 
triply ionised carbon produced. 

Finally, we briefly examine the observational conse- 
quences for studies of the IGM metallicity using C iv and 
Siiv absorption at z ~ 3. We perform a pixel optical depth 
analysis on our synthetic spectra, and find that the predicted 
UVB hardness fluctuations will have little impact on obser- 
vations compared to spatial variations in the IGM metallic- 
ity. We conclude that the lack of any observed fluctuations 
in the tsuv/tciv ratio does not provide a stringent limit on 
the non-uniformity of the UVB spectral shape, and in par- 
ticular does not preclude the possibility of He n reionisation 
completing around z ~ 2 — 3. On the other hand, we con- 
firm the observed scatter in the IGM metallicity inferred 
from C iv and Si iv absorption z ~ 2 — 3 is likely to be in- 
trinsic, reinforcing its potential as a powerful constraint on 
intergalactic metal enrichment scenarios. 
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